Geometric properties of graph layouts optimized for greedy navigation 
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The graph layouts used for complex network studies have been mainly been developed to improve visual- 
ization. If we interpret the layouts in metric spaces such as Euclidean ones, however, the embedded spatial 
information can be a valuable cue for various purposes. In this work, we focus on the navigational properties 
of spatial graphs. We use an recently user-centric navigation protocol to explore spatial layouts of complex 
networks that are optimal for navigation. These layouts are generated with a simple simulated annealing opti- 
mization technique. We compared these layouts to others targeted at better visualization. We discuss the spatial 
statistical properties of the optimized layouts for better navigability and its implication. 

PACS numbers: 89.40.-a, 89.75.Fb, 89.75.-k 



I. INTRODUCTION 

One of the quintessential properties of spatial structures is 
to support human navigation, or give cities, buildings etc. their 
navigability . Navigability in purely topological (not spatially 
embedded) graphs or networks have been studied, especially 
during the last decade of complex network research [ 1 j. Navi- 
gability in more relevant context of spatial networks |2|, how- 
ever, is a more recent topic [3). As a contribution to this 
topic, we have introduced greedy spatial navigation (GSN) as 
a simple probe of the effects of human cognitive limitation to 
navigation in non-familiar environments |@]|5]. The essence 
of GSN is the quantification of navigability for given spatial 
structures. In this work, we focus on the reverse problem of 
giving spatial structures for better navigability. 

In our previous work, we made the observation that a 
spring-embedding layout (for better visualization) actually 
helps the greedy navigators to find better routes compared 
to a random layout [4;]. The difference is occasionally very 
large, depending on the network structure. This "side effect" 
comes from the property that layouts for visualization usu- 
ally put topologically close vertex pairs geometrically close. 
In this paper, we study the layouts optimized solely for the 
better greedy navigability. In particular, we focus on geomet- 
ric characteristics of optimized layouts compared to the afore- 
mentioned visualization oriented layouts. The optimization is 
done by simulated annealing (SA) process |6). This will be 
further discussed in the next section. 



II. LAYOUT OPTIMIZATION PROCEDURE 

We begin the optimization procedure by assigning a ran- 
dom position for each vertex independently, inside square of 
unit length. The object function to minimize is the average 
steps greedy navigators take for all the vertex pairs as source- 
target (s-t) ones (denoted as d g in Ref. [5 |). We use the SA [6 1 
technique that repeatedly applies the heating and quenching 
processes. At each time step, the position of a randomly 




FIG. 1: (color online) Examples of the optimal (a) and KK layout (b) 
of the BA model. The GSN pathway is 3.85 (4.79) for the optimal 
(KK) layout, respectively. 



chosen vertex, whose coordinates are (xo,yo), is relocated to 
(xo + Ax, yo + Ay) where Ax and Ay are randomly drawn from 
the interval [-/, l\ (I » 1) uniformly. In the heating (quench- 
ing) process, such a trial movement is accepted if d g is de- 
creased, while it is accepted with the probability /?high (.Plow) 
otherwise, respectively. This is similar to successful heuristic 
methods in combinatorial optimization |7|. The trial move- 
ment in the heating process is repeated for 7# times in the 
unit of Monte Carlo (MC) steps, where one step is defined as 
N (the number of vertices) trial movement. On the other hand, 
in the quenching process, the trial movement is repeated until 
consecutive Ti times (again, in the unit of MC steps) of rejec- 
tion occurs. Overall, these consecutive heating and quench- 
ing processes as a single session are repeated Thl times in 
total. Finally, during this procedure, the layout L m ; n corre- 
sponding to the minimum d g value up to present is recorded 
and constantly updated if a new minimum d g occurs. With the 
method, after the SA procedure ends, we obtain the approxi- 
mated optimal layout for GSN up to the moment. 
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FIG. 2: (color online) A typical time series of d g in the unit of MC 
steps ?mc, in case of the BA model used in Fig.[T] along with the real 
shortest path length d, d g for the optimal layout L min [Fig.[TJa)] and 
the KK layout [Fig. |TJb)], and the moment of L min denoted as the 
vertical line. 



III. RESULTS 

To start our exploration, we set I — 10, phi g h = 0.2, 
Plow = (completely frozen phase), 7# - 500, Ti = 100, and 
Thl - 10. Due to the computational complexity we present 
the results using the graphs with rather small sizes: N = 50, 
in case of Barabasi-Albert (BA) CD, Holme-Kim (HK) 0, 
Watts-Strogatz (WS) [ 1 1 model graphs, where the average 
degree is set as k — 2 for BA and HK and k = 4 for WS 
graphs (triangle formation probability for HK model is 1 and 
the rewiring probability for WS model is 0.1). As an exam- 
ple of real-world graph, we analyze the social network of the 
oft-studied Zachary karate club with N = 34 ifTTI . As repre- 
sentative examples of completely regular structures, the two- 
dimensional square lattice (2D square) with the open bound- 
ary condition with N = 7 x 7 = 49 and the one-dimensional 
ring (ID ring) with N = 50 and k — 2 (connected only with 
the nearest neighbors) are analyzed as well. A typical exam- 
ple of L m i n and time series of d g in case of the BA model is 
illustrated in Figs.[TJa) and|2] 

In our further analysis we first investigate the angle (defined 
as the undivided angles between edge pairs attached to each 
vertex — in Fig. [Tjthis is means we include 0\ and 9i but not 
9\ + 62) distributions of optimal layouts for GSN, in compari- 
son to the spring-embedding Kamada-Kawai (KK) layout for 
the purpose of visualization |T2"I . as shown in Fig. [3] An ex- 
ample of L m ; n and time series of d g in case of the KK layout 
is shown in Figs.[TJb) and [2] Even if it is demonstrated that 
the KK layout is helpful for GSN compared to the random 
navigation [4|, the optimal layout shows significantly differ- 
ent angular profiles to the ones in KK layout. In particular, 
much sharper angles dominate the former case, in contrast to 
the latter case composed of a certain range of characteristic 
angles for better visualization. Among the optimized layouts, 
the graphs with hubs (vertices with large numbers of neigh- 
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FIG. 3: (color online) Distribution of angles in the optimized (a) and 
KK (b) layout. At least 39 graph ensembles are used to average for 
all the cases. 
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FIG. 4: (color online) Average step decreased along the GSN paths 
in terms of the relative position / along the pathways, in the opti- 
mized (a) and KK (b) layout. At least 39 graph ensembles are used 
to average for all the cases, and the horizontal axis is equipartitioned 
into 10 bins. 
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FIG. 5: (color online) An illustration of the definition of deviation 
from the straight line for a s-t pair. Each intermediate vertex's lo- 
cation in the GSN pathway of length 4 (arrows) is compared to the 
corresponding intermediate points (green dots) equally spaced in the 
straight line (red line). The deviation for each intermediate point 
is defined as the Euclidean distance between the two points (blue 
dashed lines), in the unit of the straight line. 



bors) such as the BA, HK, and Karate Club graphs, are obvi- 
ously dominated by more sharper angles attached to the hubs 
than the WS graph case. 

For a deeper investigation of GSN pathways, we analyze 
the average distance (number of edges in the shortest path to 
the target) decreased along the GSN pathway. We denote this 
quantity as (Ad(f)) and present in Fig. |4] where the relative 
position in a pathway is denoted as the fraction / = n/d g 
(n = 1, • ■ • ,d g — 1) for the n'th intermediate vertex. As seen, 
for all the graphs, \(Ad(f))\ is larger for the KK layout than 
for the optimal layout in the early stage of the GSN pathways. 
However, the optimal pathway shows larger \(Ad(f))\ soon 
afterwards, which leads to better performance (shorter GSN 
pathways). Another aspect is the peak of {Ad(f)) for both the 
KK and optimal layouts, indicating the characteristic interme- 
diate peak (point of inefficiency) in the GSN pathways. 

In contrast to the random graph models (BA, HK, and WS) 
and Karate Club graph, in case of the regular lattices (2D 
square and ID ring) the KK layout provides the exactly op- 
timal layout for GSN, where the topological square lattice 
is mapped to a geometrical square lattice with almost right 
angles. This can be seen from the peaks of Fig. [3jb). The 
topological ID ring network is laid out as a circle (with al- 
most straight angles as shown in the peak in Fig.[3jb)], and by 
(Ad(f)} = -1 regardless of / in Fig.gb)). The SA algorithm, 
however, shows a suboptimal performance for such regular 
structures as shown in Figs. [3ja) and[4ja), in contrast to the 
KK layout very good at detecting such regular structures. 

Another way to explore the GSN pathways is to measure 
their straightness, by comparing the pathways to the straight 
line connecting s-t pairs. Similar to the analysis of (Ad(f)), 
we quantify the local straightness as a function of /, as illus- 
trated in Fig. [5] The deviation from the straight line connect- 
ing a s-t pair, represented as a vector v Jf from s to t, at / is 
denoted as dev(/) — the Euclidean distance between the inter- 
mediate vertex at / and the point at the same fraction / on 
the straight line. Note that all the distance measures for each 
s-t pair are scaled with respect to the length of the straight 
line (|v. sf | = 1), to facilitate the comparison between all the 
different paths and layouts. Furthermore, each vector dev(/) 
(whose direction is from the point on the straight line to the 



vertex on the GSN pathway) for the distance dev(/) is decom- 
posed into the component parallel to the straight line devy(/) 
and the one perpendicular to it dev ± (/). Due to the intrin- 
sic directionality from the source to the target, dev||(/) can be 
positive or negative depending on the angle between \ st and 
dev(/), and we take the positive value of dev ± (/) for conve- 
nience. 

For the optimized and KK layouts applied to the graphs we 
have used, values of average deviation from the straight lines 
as a function of / are presented in Fig. [6] — (dev(/)), (devy(/)), 
and (dev ± (/)) separately. From the results, we observe that in 
most cases the GSN pathways with KK layouts show much 
less deviation from the straight lines than with the optimized 
layouts. At a first glance, it may seem counterintuitive, be- 
cause the GSN routes are shorter in the optimized except for 
2D square and ID ring cases. It is, however, reasonable if 
one recalls the fact that the path length is defined as the hop- 
ping distance (number of edges along a pathway), so that op- 
timized layouts can have much longer edges compared to the 
KK ones [compare Figs.[T|a) and (b)]. Such longer edges am- 
plify the entire length scale of optimized layouts compared to 
KK ones. Therefore, in case of hopping-distance-based opti- 
mized layouts, the absolute values of deviation themselves are 
not quite comparable to different types of layouts. However, 
for a single type of layout, we can compare the effects of dif- 
ferent graph topologies. For instance, the (purely topological) 
BA model shows the largest deviation in most cases, related 
to the amount of embedded geometric information 0. The 
circular trajectories of GSN pathways in the ID ring case are 
clearly reflected in the concave curves of dev(/) and devj_(/) 
and the slightly negative (positive) dev||(/) for small (large) /, 
respectively, for both optimized and KK layouts. It would be 
interesting if we adopt the alternative definition of path length, 
for instance, the sum of Euclidean distance along the path lfT3l 
and compare the results. 



IV. SUMMARY AND DISCUSSIONS 

We have investigated the properties of optimized spatial 
graph layouts for GSN generated by a simple SA process. 
From the observed geometrical measures, it is shown that, 
in general, the optimized layouts are characterized by sharp 
angles and point of inefficiency in the middle of the GSN pro- 
cesses. These properties are qualitatively different from the 
layouts for better visualization, namely the KK ones, show- 
ing the dominating intermediate angles within characteristics 
ranges. A closer inspection of the navigational routes also re- 
veals that the inefficiency of KK layouts soon follows its ini- 
tially better performance compared to the optimized ones. In 
other words, our simulation shows that it is possible to gener- 
ate GSN-friendly layouts other than the visualization-friendly 
layouts by just taking the simple SA optimization process. 

Our simple SA optimization process, however, also shows 
its limitation by yielding the suboptimal results for completely 
regular structures such as the ID ring and 2D square lattice, 
where the KK layouts happen to coincide with the exact solu- 
tion of optimization for GSN. The SA optimization certainly 
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FIG. 6: (color online) Average deviation from the straight line connecting s-t pairs, for the relative position x in the GSN pathways and straight 
lines. Upper (lower) panels correspond to the optimized (KK) layouts, respectively. (dev(/)) [(a) and (d)] is decomposed into (devy(/)> [(b) 
and (e)] and <dev ± (/)) [(c) and (f)] with respect to the straight line. At least 39 graph ensembles are used to average for all the cases. The 
horizontal axis is binned into 10 equidistant intervals. 



gets close to the exact optimal layout, but not perfectly, at 
least in our simulation setting. Besides such possible subop- 
timal performances, the observed properties of optimized lay- 
outs for better navigation can, we believe, give valuable hints 
for constructing various spatial structures in practice, e.g., the 
urban planning and architecture lfl4l . Adopting more sophis- 
ticated optimization processes and studies on diverse graph 
structures would be a natural candidate for the future work, 
for even better understanding of the tripartite relationship of 



topology-geometry -navigability. 
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